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ABSTRACT 

We estimate the average radio-AGN (mechanical) power deposited into the hot atmospheres of galaxy clus- 
ters over more than three quarters of the age of the Universe. Our sample was drawn from eight major X-ray 
cluster surveys, and includes 685 clusters in the redshift range 0.1 < z < 0.6 that overlap the area covered by 
the NRAO VLA Sky Survey (NVSS). The radio-AGN mechanical power was estimated by scaling the radio 
luminosity of central NVSS radio so urces using the relation between the radio synchrotron luminosities and 
X-ray cavity power measurements of Cavagnolo et al. (2010). We find only a weak correlation between radio 
luminosity and cluster X-ray luminosity across the sample. This trend is driven primarily by the most distant 
clusters, where the detection fraction and average radio powers are higher in the most luminous X-ray clusters 
at redshifts at or above z = 0.3. The average AGN mechanical power of 3 x 10 44 erg s _1 exceeds the X-ray 
luminosity of 44% of the clusters in our sample, indicating that the accumulation of radio-AGN energy is sig- 
nificant in these clusters. Integrating the AGN mechanical power to redshift z = 2.0, using simple models for 
its evolution and disregarding the hierarchical growth of clusters, we find that the AGN energy accumulated 
per particle in low luminosity X-ray clusters exceeds 1.0 keV per particle. This conservative estimate is com- 
parable to the level of energy needed to "preheat" clusters, indicating that continual outbursts from radio-AGN 
are a significant source of gas energy in hot atmospheres. Our result implies that the supermassive black holes 
in brightest cluster galaxies that generated this energy did so by accreting an average of ~ 10 9 M Q over time. 
This is comparable to the rapid level of growth expected during the quasar era, and in extreme cases may have 
lead to the development of ultramassive black holes. 

Subject headings: Galaxies: clusters: general; Galaxies: clusters: intracluster medium; Galaxies: quasars: 
general; X-rays: galaxies: clusters; Radio continuum: galaxies 



1. INTRODUCTION 

Models for the formation and evolution of cosmic structure 
generally invoke some heating mechanism to prevent catas- 
trophic gas co oling and excessive star formation in massive 
galaxies (e.g.. ISiiacki & Springel 2006). In galaxy clusters, 
the same mechanism may generate the excess entropy respon- 
sible for deviations fro m the self-similar s caling relations ex- 
pected otherwise (e.g.. lMarkevitchlfl998l) . For example, the 
Lx - T rel ation for galaxy groups is steeper (Lx oc T 2 6 ; 
Markevitchll 19981) than for clusters (L x oc r 2 : lKaiseii[T986t 
Arnaud & Evrardlll999l) . Excess entropy is also revealed by 
flatter core en tropy profiles in galaxy groups than in massive 
clusters (e.g.. IVoit & Donahue 2005). Furthermore, in cool- 
ing core clusters, the heating rate needs to be related to the 
high rate of radiative cooling that w as previously thought to 
cause cooling flows (Fab ianl ll 9941) . The heating cannot be 
too effective, since cooling co res are found in a large fraction 
of local X-ray clusters (e .g., [Mittal et al. 2009; Hudso net al.l 
120 10t ISantos et"aT1l2010l) . but it must be sufficient to explain 
the scarcity of cooling gas that would be expected t o accom- 
pany strong cooling flows (e.g. JPeterso n et al. 20031). 

One of the most promising sources of this non-gravitational 
energ y is active galaxy nuclei (AGN; e.g., McNama ra et al.l 
2000). The questions remain of when, how and and how 
much AGN energy is distributed into their environments (e.g., 
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IShort et al 112012b lYoung et al.ll201 ll) . The preheating model 
of Kaiser (1991) proposes that energy injected into the inter- 
galactic medium at high redshifts explains the observed de- 
partures from self-similar scaling relations. IWu et al.l (120001) 
found that the minimum excess energy required to break self- 
similarity is ~ 1 keV/particle. At high redshifts, many AGN 
are in the radiatively efficient "quasar" mode (ICroton et al.1 
2006), when high AGN accretion rates are promoted by the 
high galaxy merger rate. Although most of the energy out- 
put of these AGN is radiated away, quasars are so power- 
ful that only a small fraction of this energy is required to 
produce the excess entropy in hot atmospheres. By con- 
trast, there are far fewer quasars at lower redshifts, but X- 
ray observations have revealed that AGN in "radio mode" 
deposit significant amounts of energy into their hot atmo- 
spheres. The total energy output of AGN in radio mode is 
generally less than in quasar mode. Nevertheless, a large 
proportion of this emerges as mechanical energy in jets and 
simulations suggest that the radio mode heating is neces- 
sary, in addition to the preheating, to suppress c ooling flows 
in clusters and star formations in galaxies (e. g., iBower et al.l 
120061: ICroton et al.l 120061: ISiiacki et al.l 120071) . Radio AGN 
hosted by cluster central galaxies in the local Universe have 
been shown to deposit enough power to prevent rapid cool- 
ing and star formation in the centers of many clusters (e.g., 
iBirzanet al.l 12004 iBest et all 120071) . Supported by the cor- 
relation between radiative cooling rate s in clusters and the 
radio power of the c entral AGN (e.g.. iRaffertv et al] 120061; 
iDunn & Fabianl I2006h . the power output of the AGN is be- 
lieved to be coup led to the cooling rate of th e hot gas in 
a feedback loop dMcNamara & NulseiJ 120071) . Neverthe- 
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less, some powerf ul AGN apparen tly reside in non-cooling 
core clusters (e.g.. ISun et al.ll2007l) . Although these systems 
lack large scale cooling flows, accretion from small coronae 
around their AGN ca n support powerful radio sources (e.g., 
IHardcastle et al.l 120071) . Radio AGN in the non-cooling core 
clusters have been shown to contribute si gnificantly to the ex- 
cess entropy of less massive clusters (e.g., Giodini et al .l2010t 
iMa et al.ll20ll . 

In this paper, we focus on estimating the average mechan- 
ical power output of radio AGN in clusters out to z ~ 0.6, 
corresponding to a look-back time of about 5.7 Gyr. At 
these redshifts, it is difficult to estimate directly from X- 
ray images the amount of energy deposited by AGN in the 
intracluster medium (ICM). The only systematic search for 
X- ray cavities at these redshifts w as undertaken recently 
by Hlavacek-Larrondo et al. (2012a). They concentrated on 
identifing cavities in the most luminous, and bright cooling 
core clusters and so could not quantify the average AGN 
output of clusters overall. The question of how much en- 
ergy AGN contribute to the ICM at higher redshifts is sig- 
nificant because AGN activity increa ses with redshift (e.g., 
iMartini et aT1l2009l) . lHart et all (1201 11) suggest that the power 
injected into clusters by radio AGN at redshift 1.2 is substan- 
tial, a factor of 10 greater than injected locally. In addition, 
the fraction of cooling core clusters appears to evolve with 
redshift, so that man y fewer large cooling cores are found be- 
yond z ~ 0.3 (e.g. 
but see iSantos et al.l 



Santo s et al.|[2Q10t ISamuele et al]|20lU 
2012J). If jet power is coupled to cool- 



ing power by a feedback loop, this suggests that the mean jet 
pow er in high - redshift clusters should be reduced. 

In lMa et alldMTllMlNTll hereafter), we estimated the av- 
erage mechanical energy deposited by radio AGN in galaxy 
clusters using the clusters in the 4 00 Square Deg ree Cluster 
Survey (400SD, 0.1 < z < 0.6: iBurenin eTal1l2007l) and 
the radio sources i n the NRAO VLA Sky Survey (NVSS; 
ICondon et al.lll998l) to show that the AGN heating in radio 
mode could also contribute significantly to the energy budget 
of clusters and groups. We found that 30% of the clusters 
showed radio emission within a projected radius of 250 kpc 
and above a flux threshold of 3 mJy, despite t he declining 
numb e rs of cooling core c lusters in the 400SD (Santo s et al. 
"20TO [ ISamuele et all |20TTI see also McNamara & Nulsen 
2012t iMann & E beling 2012). The average jet power of the 



central radio AGN is approximately 2 x 10 ergs . As- 
suming that the current AGN input power remains constant to 
redshifts of 2, the energy input per particle would be at least 
0.4 keV within Psoo- In addition, we found no significant 
correlation between the radio power, i.e., the mechanical jet 
power, and the X-ray luminosities of clusters in the redshift 
range 0.1 - 0.6. This implies that the mechanical jet power 
per particle is higher in clusters with lower masses. However, 
within this single flux-limited cluster survey, the X-ray lumi- 
nous clusters are also the clusters with the highest redshifts. 
Thus, we could not distinguish redshift evolution from lumi- 
nosity dependence for AGN heating. In the present study, we 
try to break this degeneracy using a composite sample from 
eight X-ray cluster surveys. 

The method used to estimate jet mechanical powers from 
the radio powers of cluster central galaxies is reviewed in Sj2] 
The composite cluster sample is introduced in $3] In SjU we 
examine the correlation between the power of a radio galaxy 
and the X-ray luminosity of its host cluster. Ogives estimates 
for fractions of clusters with a central NVSS source and av- 
erage radio powers as functions of redshift and X-ray lumi- 



nosity. The evolution of cluster radio power is discussed in 
Sj6] The energy per particle deposited in groups and clusters 
since redshift 2 is estimated in $7] $8] contains some discus- 
sion of the calculation of average AGN jet power and £|9] is 
the summary. We adopt a ACDM cosmology with ho = 0.7, 
tt A = 0.7, and Vt m = 0.3. 

2. AGN HEATING: MECHANICAL POWERS OF JETS 

X-ray cavities provide clear evidence of the interaction be- 
tween AGN jets and the hot atmospheres of clusters. Power 
from AGN can be distri buted into the ICM through sev - 
eral channels (reviewed in McNamara & Nulsen 2007, 2012), 
e.g., shock fronts dNulsen et al.l l2005allbl) and sound waves 
(Fabian et al. 2006) driven by AGN jets. The minimum en- 
ergy required to create a cavity can be estimated using simple 
assumptions. The enthalpy, H, of a cavity is equal to the sum 
of its thermal energy and the work required to excavate it un- 
der constant pressure. If the plasma filling the cavity is pre- 
dominantly relativistic, this is H = 4pV, where p is the pres- 
sure and V is the volume of the cavity. Assuming a time scale, 
r, to inflate the cavity, the mean power required to inflate the 
cavity is at least Pj et ~ ApV/r, which provides an estimate of 
the jet power. The time scale, r, is commonly estimated us- 
ing the terminal velocity of the buoyantly ris ing bubbles (e.g., 
IBirzan et al.ll200l 12001 iDunn et al.l 1200 51). Such measure- 
ments of the jet power require deep, high resolution X-ray 
data to determine the volume and pressure of the cavities, so 
that the measurements cannot currently be conducted for a 
large statistical sample. Nevertheless, a correlation between 
the radio p ower of the AGN and the jet power was demon- 
strated by IBirzan et al.l d2"0 04) and i mproved bv IBirzan et al.1 
(2008). ICava"gnolo et al l d2010l) . and |Q' Sullivan et al.l d201 C 
Using this correlation, we can estimate the minimum power 
necessary to inflate radio lobes from the radio power of the 
central AGN. This procedure provides a practical means to 
estimate the energy deposited by radio AGN in a sample large 
enough for a statistically meaningful analysis. 

In this work, w e estimated jet powers u sing the Pi .4 - Pj et 
scaling relation of Cava gnolo et al.1 (120101) . 

log P jct =0.75(±0.14)logPi. 4 + 1.91(±0.18), (1) 

where Pj 0t is in units of 10 42 ergs~ 1 and Pi ,4 is the ra- 
dio power at 1.4 GHz in units of 10 40 ergs _1 . The scat- 
ter in the correlation between jet power and radio luminos- 
ity is (T1.4 = 0.78 dex. Although measurement errors con- 
tribute to this scatter, it is dominated by intr insic variations 
in radio source properties (IBirzan et al.ll2008l) . The relation- 
ship in Equation is determined over 7 decades in Pi 4 
(10 37 - 10 44 ergs _1 ), for systems ranging from the nuclear 
radio sources of Brightest Cluster Galaxies (BCGs) in cool- 
ing core clusters with Lx up to 10 45 erg s _1 to the low -power 
radio sources in galaxy groups, with Lx of approximately 
lO^ergs " 1 . Note that there are only three sources in the 
sample of Rlavagnolo etail (120101) with P a > lO^ergs" 1 
and the relation in Equation (fTJ may overestimate Pj 0t for 
them. This is discussed further in £13.21 

In contrast to the cavity powers, "beam" powers of radio 
AGN have been estimated based on ly on radio data (e.g., 
lO'Dea et al.l l2009t iDaly et al.1 120TI lAntognini etail I2012L 
and the references therein). Cavity powers and "beam" pow- 
ers provide largely complementary means to estimate the 
AGN jet power, since cavities are mostly associated with FR I 
radio sources, whereas the beam powers are only measured 
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for FRII ra dio sourc es (Fanaroff & Rilevll 19741). The tw o re- 
cent papers lDalv et al.l (12012b and lAnto gnini etal.ld2012l) dis- 
cuss the relationship between beam power and radio power for 
FRII sources. For a radio source with Pi. 4 ~ 10 40 ergs -1 , 
their beam powers are an order of magnitude greater than the 
cavity power calculated using Equation (TJi. In addition, the 
slopes they find for the rel ationship between beam power and 
radio power (0.84±0 14in lDalv et al.l2012l and 0.95±0.03in 
Antog nTni et al.ll20l2l) are steeper than given by Equation ((TJ. 

Nevertheless, because of the large scatter in these rela- 
tions, the slopes diffe r from the relation Equation Q] from 
(ICavagnolo et al.l2010l) by less than 2a. It is beyond the scope 
of this paper to resolve these inconsistencies. Here, we use the 
relation of Equation ((TJ, since it provides a more conservative 
jet power. If the jet powers, particularly those for powerful 
radio sources, are underestimated as suggested by the beam 
power - radio power relations, our estimates of average jet 
power are also underestimated, so that the mechanical energy 
deposited by radio AGN is even more important. 
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FIG. 1. — Bolometric X-ray luminosity versus redshift for clusters in the 
eight cluster surveys used here. Bolometric X-ray luminosities are derived 
from 0.5 - 2.0 keV ROSAT PSPC luminosities, assuming the L x - T re- 
lation of Markevitch ( 1998). Symbols for the cluster surveys are given in 
the legend. Red symbols denote clusters having an NVSS source projected 
within 250 kpc of the cluster center. The gray boxes mark the cells used to 
calculate average jet powers in S|7] 



3. THE SAMPLE 

To extend the work of IMMN1 ll we have com- 
bined e ight major X-ray c luster surveys (Figure [TJ: the 
400SD (Il3menin etal.ll2007l). the 160 deg 2 Survey (160SD 
iVikhhnin et al.l Il998t iMullis et al.l 120031), the Wide An- 
gle ROSAT Pointe d Survey (WARPS: iScharf et all [l997t 



iHorne r et al. 2008), the MAssive Cluster Surve y (MACS 
lEbeling et alj|200ll 120071 l20lot iMann & EbelT^l2012 fl the 
Brightest Cluster Survey (BCS: lEbeling et al.l 119961 2000) 
the C l usters in the Zone of Avoidance (CIZA; Ebeli ng et al 
120021 lKocevskietaTl |2007| ), the Northern ROSA T All-Sky 
Cluster Survey (NORAS; IBohringer etall 120001) , and the 
ROSAT-ESO Flux Lim ited X-Ray Cluster Survey (REFLEX; 
IBohringer et aLll200Tb . 

4 We only used the MACS subsample publi s hed in the three papers of 
lEbeling etal] ( 2007. 2010) and Mann & Ebe Knel UM% . 
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FIG. 2. — Correlation between cavity power, P C av, estimated as ApV/r, 
and radio power, Pi. 4. Clusters in our sample with known X-ray cav- 
ities are plotted as filled symbols. P C av values for the m are taken 
from Hlavacek-Larrondo et al. (2012a, filled green circles) and Birzan et al. 
120081 filled black triangles). "Possible" cavities from the sample of 
Hlavacek-Larrondo et al. (2012a) are marked with a cross to indicate that they 
are less secure detections. Open triangles and circles d enote other clusters 
used by Birzan et al. (2008) and Cavagnolo et al. 12010) respectively. The 
solid line is the scaling relati on o f Equation CD- The dotted line shows the 
saturation level discussed in i|3.2l 



The first three surveys were compiled from serendipi- 
tously detected clusters in targeted ROSAT PSPC observa- 
tions, while the other five cluster survey s are based on th e 
ROSAT All Sky Survey catalogue (RASS: fVbges et alJH999h . 
Many clusters are recorded in more than one of these sur- 
veys. These were identified by having centroid offsets of 
less than 2'. Most overlapping identifications are between 
the three serendipitous surveys, or NORAS and BCS. Of the 
223 160SD clusters, 101 are included in 400SD, 40 of the 
141 WARPS clusters are included in 400SD, and 213 of the 
484 NORAS clusters are included in BCS. Redshifts for the 
overlapping entries are mostly consistent within 10% between 
catalogs. For the overlapping entries, we used the redshift, 
centroid, and X-ray luminosity from the highest priority sur- 
vey. Priority order is as listed abovfl i.e. 400SD, 160SD, 
WARPS, MACS, BCS, CIZA, NORAS, and REFLEX. For 
the overlapping entries, the multiple redshift measurements of 
18 clusters are inconsistent, with |Az/z| > 10%. These clus- 
ters were excluded from our sample to ensure the accuracy of 
our redshifts and cluster identifications, although some mea- 
surement inconsistencies see m to have been reso lved in the 
literature (see the summary in Piffaretti et al. 20TJJ). 

To combine the X-ray luminosity measurements from dif- 
ferent surveys, we recalculated b olometric X-ray l uminosities 
using the L x - T relation from lMarkevitchl (119981) . The X- 
ray luminosities estimated from different entries for the same 
cluster were compared to examine the consistency of flux 
measurements in different surveys. The average difference 
in luminosity is about 10%. X-ray luminosities are used to 
estimate the number of particles in the clusters in fJ7]and so to 

5 Of the three serendipitous surveys, the most up to date, 400SD, is ranked 
ahead of 160SD, while WARPS is ranked last for its smaller sample size 
and sp arser information in the public catalog iScharf et al. 1997; Horner et al. 
12003) . Of the all-sky surveys, the rankings of MACS, CIZA, and REFLEX 
are insignificant because they have so few overlaps with the other surveys. 
For the two remaining surveys, BCS and NORAS, the redshifts and X-ray 
luminosities are comparably reliable. We ranked BCS higher because of our 
familiarity with the BCS work. 
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TABLE 1 

Fraction of clusters with NVSS sources" 



Fig. 3. — Left: NVSS radio flux for AGN projected within 250 kpc of a cluster center versus 0.5 - 2keV X-ray flux of the cluster. Right: R adio power at 1 .4 
GHz f or the same sources versus bolometric X-ray luminosity of the host. Clusters identifi ed as having stro ng cooling cores (t coo \ < 1 Gyr) in Cavag nolo et alj 
12009) are plotted in red. Most of the remaining clusters are not from the sample of Cavagnolo et al. (2009), so their cooling times are unknown. 

derive the average AGN energy injected per particle. Because 
uncertainties in Pj 0t due to the scatter in the scaling relation 
of Equation ([T|i dominate in estimates of the average radio jet 
power, uncertainties in the X-ray luminosities are not an issue 
in this study. 

Collectively, these eight X-ray cluster surveys contain 1032 
clusters in the declination range —38° < S < 68°, an area that 
is well covered by the NVSS. The background source densitjQ 
for the sample is pbkg > 35deg -2 . We focused on clusters 
in the redshift range 0.1 to 0.6. The upper redshift limit is 
set by limited sampling and the lower limit is set to avoid the 
complexity of radio flux measurements for resolved sources. 
After the redshift cuts, 685 clusters remain in the range of 
bolometric X-ray luminosities 3 x 10 43 to 15 x 10 45 ergs -1 . 



Survey 




0.1 < z < 0.3 


0.3 < z < 0.6 




Ncl 


f H 


iR.hif 


N c i 


f b 
iR.hif 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


400SD 


97 


0.32 ±0.06 


0.09 ± 0.03 


53 


0.25 ±0.1 


160SD 


80 


0.32 ±0.07 


0.21 ±0.06 


63 


0.29 ±0.1 


WARPS 


55 


0.33 ± 0.09 


0.16 ±0.06 


47 


0.32 ±0.1 


BCS 


131 


0.55 ±0.08 


0.28 ± 0.05 


9 


0.44 ±0.3 


MACS 









65 


0.48 ±0.1 


CIZA 


35 


0.63 ±0.17 


0.29 ±0.10 
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NORAS 


191 


0.55 ±0.07 


0.26 ±0.04 


23 


0.13 ±0.1 


REFLEX 


122 


0.59 ±0.09 


0.31 ±0.06 


7 


0.86 ±0.5 



3.1. Radio Sources in Clusters 

Following the analysis in M MN1 ll we cross-matched the 
coordinates of the clusters with radio sources in the NVSS 
catalogue. For our sample of 685 clusters, 357 have NVSS 
radio sourc es above a flux li mit of 3 mJy proje c ted within 
250kpc. IMMN1 11 (see also iLin & Mohrl l2007t iBest et alj 
2007) showed that the density of radio sources at the center of 
the clusters is much higher (~ 2Mpc -2 ) than at larger radii 
(0.3 Mpc -2 ), so the probability that these central sources are 
not associated with the clusters is small. The total expected 
number of background contaminated clusters is 25, i.e., 7% 
of the 357 clusters. There is little to be gained from using a 
smaller aperture due to the large uncertaintie s in the c luster 
coordinates determined from ROSAT data (see MMN1 1 for a 
more detailed discussion based on 400SD clusters). 

In Table [T| we show the number of clusters in each survey 
and the fraction of clusters with radio sources for two red- 
shift ranges, 0.1 < z < 0.3, and 0.3 < z < 0.6. Here, the 
cluster samples are defined by the same criteria, i.e., declina- 
tion range and NVSS background density, as the composite 
sample (see Sj3]), but redundant entries for a cluster in the dif- 
ferent surveys are not excluded. For the lower redshift range 
(column 3), the radio source fractions for the three serendip- 
itous surveys, at ~ 32%, are consistently lower than those 
for the five all-sky surveys, which all exceed 55%. For the 
higher redshift range (column 6), the situation is similar, al- 
though the variations between the all-sky surveys are greater 
due to large statistical errors. These differences between the 

6 Background source density, Pbkg> is measured in an annulus extending 
from 2 to 5 arcmin from each cluster, with a flux limit of 3mjy. 



a Uncertainties in the fractions are calculated assuming Poisson errors in the 
counts. 

k fR.hif is the fraction of clusters having a central radio source with Pi. 4,11m > 
3.8 X 10 40 erg s -1 , the radio power of a 2 mJy source at z — 0.6. 

serendipitous and the all-sky surveys are probably due to the 
correlation between radio source fraction and the X-ray lumi- 
nosity of a host cluster, since the clusters in the all-sky surveys 
are generally more luminous than those in the serendipitous 
surveys at a given redshift. This correlation is discussed in 
S|5] In order to provide a fair comparison of radio source frac- 
tions at different redshifts, fR.hif in column (4) is the frac- 
tion of clusters with a central radio source more powerful 
than Pi. 4 = 3.8 x 10 40 ergs _1 , the power cut defined for 
the high redshift sample. Column (6) gives the same fraction 
for the higher redshift sample, showing that these fractions 
are marginally greater for the higher redshift range. 

3.2. Cavity Powers 

Some o f our sample clusters were show n to have X-ray 
cavitie s bv lHlavacek-Larrondo et al] d2012al) and lBirzan et alj 
(2008). Figure [2] shows estimates of jet power, P cav = 
4pV/rbuoy> for these, plotted against their radio powers, to- 
gether with the s caling relation of E quati on CD), and the 
remain ing data of Birzan et al. (2008) and Cavagnol o et alj 
(2010) used to establish the scaling relation. The new data 
(filled green circles) are consistent with the scaling relation 
and their scatter is similar to that for the data used to establish 
the scaling relation. 

The scaling relation (solid line) may be overestimating Pj ct 
for the five powe rful radio sources wit h Pi . 4 > 10 42 ergs~ 1 i n 
Figure|2](but see lDalv et aT]|2012l an dlAntognini et alj|2012[) . 
A similar result is seen in Figure 1 of lCavagnolo et al] d2010l) . 
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This departure could arise if a significant fraction of their 
radio synchrotron power is generated by "hot spots" which 
are absent from the FR I radio sources that form most of 
the scaling relation (e .g., 3C295: lHarris et al] 120001 Cygnus 
A; IWilson et al.ll2006T) . This synchrotron flux should be ex- 
clude d before applying the scaling relation. In some cases 
(e.g., lFormanetalJl2005t lLal et all I20T0I) shock fronts may 
also contribute significantly to the power deposited into hot 
atmospheres. In principle, the energy associated with shock 
fronts should also be taken into account when the Pj Ct is mea- 
sured. However, high-quality X-ray data and careful data 
analysis are necessary to detect shock fronts and to estimate 
the associated energy, which would be a major undertaking 
for a large data set. For consistency, values of P cav in Fig- 
ure |2] were estimated using only the enthalpy and the cavity 
inflation time. 

The number of powerful sources in Figure|2]is too small to 
place a strong constraint on the slope of the scaling relation at 
the high end. So it is unclear whether and to what degree we 
may be overestimating the jet power in these sources. How- 
ever, this issue is crucial for to making reliable estimates of 
-Pjet in If the scaling holds, the mean power output of 
short-lived but powerful radio sources could rival the level of 
normal radio-AGN heating over time. We therefore take two 
approaches to calculating the mean power. First, we excluded 
the most powerful radio sources, so that the resulting average 
cavity power places a lower limit on the true average cavity 
power. Second, we set Pj Ct for the powerful sources, assum- 
ing that the jet power saturates at the constant value of Pi .4 = 
2.4 x 10 41 ergs -1 for high radio powers, i.e, Pj C t(Pi.4 > 
2.4 x 10 41 ergs -1 ) = 1.07 x 10 45 ergs -1 (dashed line in 
Figure I2J- The saturation level is set to the mean value of 
log P C av for the 5 most powerful radio sources in Figure [2] 
The two approaches are compared in ||7] 

3.3. Cooling Times 

Central cooling times for 110 clusters in our sample were 
estimated using archival Chandra data for the "Archive of 
Chandra Cluster Ent ropy Profile Tables" project (ACCE PT; 
ICavagnoloet"aT]l2009h . Briefly, ICavagnolo etai] (120091) fit 
annular spectra for each cluster to determine the cooling time 
as a function of the radius, assuming a profile for the cooling 
time of the form 



t CO oi(r) =t c0 + t 



100 



100 kpc 



(2) 



where t C Q and iioo are constants. The value for the central 
cooling time, t c g, is the cooling time t coo i used in this paper. 

4. CORRELATION BETWEEN RADIO POWER AND X-RAY 
LUMINOSITY 

The dependence of the power of a radio source on the X- 
ray luminosi ty of its hosting cluster is shown in Figure [3] 
As found by MMN 1 1 , the correlation is weak. The Kendall 
correlation coefficient for the fluxes is small, at r = 0.16, 
although it differs from zero at high significance, the prob- 
ability of getting a value this large by chance being only 
1.1 x 10~ 4 . The slope of the relationship between the fluxes 
is d log -Fnvss /d log Fx 0520 = 0.28 ± 0.05 and between the 
powers it is dlogPi^/dlogLx.boi = 0.33 ± 0.05. The sur- 
prising agreement between these slopes requires a weak cor- 
relation between distance and luminosity for the sample (Fig- 
ureQ]). Since cluster X-ray luminosity increases with mass, it 



follows that the radio power of a central AGN in clusters is 
weakly dependent on the cluster mass. 

Confining attention to clusters with cooling times shorter 
than 1 Gyr, plotted in red in Figure [3] the Kendall corre- 
lation coefficient for the fluxes is r = 0.3, with a prob- 
ability of 0.04, still significant at the 95% level. For 
these clusters, d log Pnvss/<^ log Fx 0520 = 0.47 ± 0.50 
and dlogPi^/dlogLx.boi = 0.51 ± 0.30. While there is 
marginal evidence that the radio power of a central AGN is 
more strongly dependent on the X-ray luminosity in these 
clusters, the correlation is w eaker than expected from the 
work of Raffer ty et alJ (|2006), who found a relationship be- 
tween jet power and coolin g rate in clusters. Sev eral factors 
may be at work here. First, iRaffertv et alJ ([2)06]) use cavity 
powers determined from X-ray data, a fairly diect measure, to 
estimate jet powers. Using Equation (HJ to connect radio pow- 
ers to jet powers injects significant extra scatter into the re- 
lationship between cavity power and cooling power. Second, 
the dynamic range of X-ray luminosities in Figure[3]is signifi- 
cantly smaller than that of the cooling powers in Raffertv et al. 
(2006) , tending to bu ry any correlation in the scatter. Lastly, 
Raffe rty et alJ J2006) relate the cavity power to the power ra- 
diated within the cooling radius. If feedback is at work, the jet 
power should only depend on cooling in this region. Although 
X-ray emission from within the cooling radius can be an ap- 
preciable fraction of the total X-ray luminosity of a cluster, 
the correlation is diluted by X-ray emission from larger radii. 

4.1. Distribution of Radio Powers for Cluster Central AGN 

If AGN feedback prevents cooling and star formation in 
cluster central galaxies, then a high cooling rate implies a 
high AGN power; thus, the radio powers of strong cooling 
core clusters are expected to be greater. 

In Figure HI we compare distributions of radio power for 
the central AGN to examine whether the clusters with cool- 
ing cores or cavities do have greater radio powers. The red 
histogram in the upper left panel shows cooli ng core clusters, 
with t coo \ < 1 Gyr from the ACCEPT data (Cavagnol o et alJ 
2009), while the blue histogram on the upper right shows non- 
cooling core clusters, with i coo i > 3 Gyr from the ACCEPT 
data. For comparison, the distribution for the entire cluster 
sample, selected as discussed in !j3] is plotted in gray. The 
distribution of radio powers in non-cooling core clusters is 
broad, including some powerful radio sources, although the 
number of clusters with t coo \ > 3 Gyr common to our sam- 
ple and the ACCEPT sample is too sma ll to provide a robus t 
resu lt. Further com plicating matters, Raff ertv et al.l (12008b 
and ICavagnolo et al.l (120081) found that radio AGN and star 
formation activity at cluster centers associated with cooling 
flows are triggered when the central cooling time falls below 
a threshold of t coo \ < 0.5 Gyr. Our X-ray data are generally 
unable to detect such a threshold, making the distinction be- 
tween cooling-flows and non-cooling flows rather uncertain. 

For a larger sample of non- or weak-cooling core clusters, 
radio powers for400SD and 160SD clusters with z > 0.3 are 
plotted in the cyan histogram of the lower right panel. The 
400SD clusters at high redshifts were fo und to be dominated 
by non- or weak-cooling core clusters by S antos et al.l ( 120101) 
and lSamuele et all ( 120111) . Clusters in the 160SD and 400SD 
cluster samples should be similar because these two surveys 
used much the same selection criteria. To avoid bias due to the 
higher cut in radio power at higher redshifts, the gray shaded 
histogram shows all the clusters of our sample with z > 0.3 
for comparison. The means of the log of the radio power for 
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these two samples are log Pi. 4 = 40.98 ± 0.08 for the 400SD 
and 160SDclusters andlogPi. 4 = 41.07±0.06 for the whole 
sample, which are consistent with one another, showing no 
evidence of an offset between the two distributions. 

Despite this, the most powerful radio sources in our sam- 
ple do tend to be associated with cooling cores. The frac- 
tion of cooling core clusters with radio sources having Pi 4 > 
10 42 ergs^Ms 13 ± 5% (6/45), compared to 10 ± 7% (2/19) 
for the non-cooling core clusters of the upper right panel and 
3 ± 2% (1/36) for the 160SD and 400SD samples in the lower 
right panel. In summary, radio sources in the cooling core 
clusters are generally as powerful as those in the non-cooling 
cores, apart from the most powerful radio sources. On the 
other hand, we have reliable estimates of the cooling rate for 
only a small sub-sample, and better coverage of deep and high 
resolution X-ray data are required for a more robust conclu- 
sion. 

The green histogram in the lower left panel of Figure |4] 
shows radio powers for c lusters identified with cavities by 
iHlavacek-Larrondo et alj d201 2ari with clusters from our 
sample at z > 0.3 for comparison. The clusters with cavi- 
ties do have slightly greater radio powers (mean log Pi. 4 — 
41.55 ± 0.22) than clusters in our sample (mean log Pi. 4 = 
41.07 ± 0.06). Furthermore, the distribution of radio powers 
for the clusters with cavities has a longer tail at high powers. 
The fraction of these clusters having a powerful radio source 
(P1.4 > 10 42 ergs" 1 ) is 30± 13% (6/20), compared to 8±3% 
(8/96) for our sample. 
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FIG. 4. — Distributions of radio power. The panels show histograms of 
Pi. 4 for various subsamples in different colors: red in the upper left panel 
for the cooling core clusters of Figure [5] blue in the upper right panel for 
non-cooling core clusters, i.e., those with cooling times greater than 3Gyr, 
cyan in the lower right panel for 400SD and 160SD clusters at z > 0.3, 
and green in the lower left pane l for the clusters with cavities identified by 
Hlavacek-Larrondo et al. 1 2012a). The gray histogram in the two upper pan- 
els shows our entire sample, as described in fTJ] and the filled gray histogram 
in the lower two panels show clusters in our sample with z > 0.3. The col- 
ored histograms are weighted by a factor noted in each panel to make com- 
parisons easier. Biweight means and their 95% confidence ranges are given 
in the same colors as the corresponding samples. 



5. FRACTION OF CLUSTERS WITH RADIO AGN 

7 Clusters with "clear " and "potential" cavities from 
IHlavace k-Larrondo et al. 12012a) are included. Three of their "poten- 
tial" cavities show no detected radio source. 
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FIG. 5. — Fraction of clusters having at least one NVSS source projected 
within 250 kpc of the cluster center. The NVSS sources are selected to have 
Pi. 4 > 3.8 X 10 40 ergs — 1 . The redshift and X-ray luminosity bins are 
defined in FigurefT] Colors are used to distinguish bins dominated by clusters 
from the all-sky surveys (green) from those from the serendipitous surveys 
(blue). Error bars show uncertainties calculated assuming Poisson statistics. 
The fraction for each bin is noted in the bin. 

IMMN1 ll found that the probability of a more luminous and 
higher-redshift 400SD cluster hosting an AGN is marginally 
higher than that for a less luminous and lower-redshift 400SD 
cluster, at the la level. In their relatively small sample, red- 
shifts and X-ray luminosities are coupled, so that the more 
luminous clusters also have higher redshifts. With our larger 
composite sample, we can examine separately how the frac- 
tion of clusters matched with NVSS radio sources depends 
on the X-ray luminosity and redshift, as shown in Figure [5] 
For computing radio source fractions here, radio sources are 
defined as having powers, Pi. 4 > 3.8 x 10 40 ergs -1 , corre- 
sponding to a radio flux of 2 mJy for a source at z = 0.6. The 
fraction of clusters having a radio source is computed for each 
bin defined in Figure Q] From the figure, the fraction gener- 
ally increases with both the redshift and the X-ray luminosity 
of a cluster. 

A significant concern is that trends in the radio fraction can 
be masked by differences between the serendipitous and all- 
sky surveys. For example, for the four redshift bins bounded 
by z = 0.14, 0.17, 0.23, 0.35, and 0.6, in the luminosity range 
3 < Lx/10 44 ergs _1 < 10, the fractions of radio sources in 
the two lower redshift bins, 0.14 < z < 0.23, which are 
dominated by clusters from the all-sky surveys, are higher 
than the fractions for the two higher redshift bins, which are 
dominated by clusters from the serendipitous surveys. This 
concern is related to the persistent question of whether the 
serendipitous surveys preferentially select non-cooling core 
clusters , while the all-sky surveys favor cooli ng core clus- 
ters (cf. lVikhlinin et alJ2006tlEckert et al.l20TTh . Under AGN 
feedback models, coolin g core clusters are more likely to host 
centra l radio AGN (e.g., Raff ertv et alj2006tlCavagnolo et alj 
2008). However, even if the serendipitous and all-sky surveys 
differ, Figure[5]suggests that this does not mask an increasing 
trend in the radio fraction with redshift and luminosity among 
the higher-redshift and more luminous clusters. This trend 
can be seen separately in the blue and green points that are 
dominated by serendipitous and all-sky surveys, respectively. 

Figure|6]shows average radio powers per cluster for the bins 
of Figure Q] As for Figure [5] a threshold on the radio power 
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of Pi. 4 > 3.8 x 10 40 ergs -1 was used to avoid spurious red- 
shift dependence in the results. The three most powerful radio 
sources in our sample, 3C 295, Hercules A, and 3C 288, are 
excluded from the averages, since they are so dominant that 
they would obscure the underlying trends. Thus, the average 
radio powers in Figure [6] are lower limits. 




FIG. 6. — Average radio power for sources projected within 250 kpc of clus- 
ter centers. Radio sources are selected from the NVSS, above a radio power 
ofPi.4 > 3.8 X lO^crgs" 1 for all redshifts. The average radio power is 
recorded in each bin. Colors have the same meaning as in Figure|3] 



6. EVOLUTION OF NUMBER AND POWER OF CLUSTER CENTRAL 
RADIO SOURCES 

The distribution of the number of radio sources per cluster, 
per unit log Pi. 4 is 



1 dN SIC (>P 1A ) 



d log Pi 



(3) 



where N aTC (> Pi. 4) is the number of radio sources with pow- 
ers greater than Pi, 4 in the cluster population of interest and 
N c \ is the number of clusters in the population. The ex- 
pected number of background radio sources for each cluster 
is subtracted from N STC and there may be more than one radio 
source in a cluster. Note that, since the normalization of ip 
gives the mean number of radio sources per cluster it may be 
greater than unity. The distribution ^(Pi.4) is plotted in the 
left panel of Figure [7] for clusters in the four redshift ranges, 
0.05 - 0.1 - 0.2 - 0.4 - 0.6. Following g] only clusters with 
luminosities in the range 3 x 10 43 < Lx < 15 x 10 45 erg s _1 
are included. Values for (^(Pi^) are shown only for pow- 
ers above a threshold corresponding to the flux limit of 3 mJy 
for the redshifts z = 0.1,0.2,0.4, and 0.6. The distribution, 
^(Pi.4), increases with redshift, evolving more significantly 
at th e higher power end, consistent with previous findings 
(e.g. JGalametz et al.ll2009trHart et al.11201 11) . 

Closely related to (p(Pia), we define </>(Pj Ct ), the number 
of radio jets per cluster per unit log Pj C t , by replacing the radio 
power Pi. 4 in Equation (0 with Pj Ct calculated from Equa- 
tion (0}, 



than Pffi is then 

$ 2 (P j 1 r)=/ 0(P jet )P j etdlogP je t. (5) 

This is plotted as a function of P^™ in four redshift ranges in 
the right panel of Figure [7] Here and earlier in Equation (O, 
d log P is calculated as the Voronoi interval for the power of 
each radio source, although, for ip( K P 1A ) 7 the data are binned 
in Pi. 4. Values are plotted only for Pj 1 ™ above a lower limit 
corresponding to the radio flux limit of 3 mJy for the different 
redshift ranges. 

6.1. Correcting for the Radio Flux Limit 
A fair comparison of the average P ]ct per cluster for differ- 
ent redshifts requires that we use the same value of P^t* 1 • For 
the whole sample, that would limit us to using the jet power 
corresponding to the radio power cutoff for z = 0.6. This 
is very restrictive and it would mean discounting the power 
input of many less powerful radio sources seen at lower red- 
shifts. Alternatively, we can estimate the average jet power 
for smaller values of Pj ™ by applying a correction factor w, 
calculated from the form of <I> Z at lower redshifts, on the as- 
sumption that the shape of $ z (Pj'™) does not evolve with z. 
This gives the correction factor 



w(z) 



<W^T(fo)] 
(*)]' 



jet 

^ZQ L- 1 jet 



(6) 



where P$£(z) is the lower limit on the jet power for redshift 
z, which is obtained by inserting the radio power correspond- 
ing to the flux limit of 3 mJy at redshift z into Equation (Q]). 
The cumulative jet power, <5> Zo , used for reference here is that 
for the redshift range of 0.05 < z < 0.1. For example, for 
the redshift bin 0.2 < z < 0.4 the correction factor is w(z ~ 
0.3) = 1.26, calculated for P${z) = 8.5 x lO^ergs" 1 and 

Pj 1 ™(zo) = 8.6 x 10 42 ergs . The correction factor, w(z), 
is used in the calculation of the average jet powers below. 



7. AVERAGE JET POWER 

The average jet powers, (Pj e t), shown in the upper panel 
of Figure [8] were estimated using a Monte Carlo method that 
accounts for the uncertainties in the radio fluxes and the pa- 
rameters in Equation ([T), the distribution of radio spectral in- 
dices, and the large intrinsic scatter (01.4) in the relation of 
Equation (fl}. Note that the log of the arithmetic mean of Pj Ct 
for a lognormal distribution is greater than the "mean" of its 
log that is given by Equation Q). For each range of Lx, the 
redshift bins of Figure [T] are chosen to distribute the clusters 
evenly between the bins. First, (Pj e t) is calculated for each 
bin in Figure Q] then this is integrated over time for a given 
range of Lx to give the time averaged mean jet power 



(P 



jet/ mt 



Z t [(P)ct(zj,Lx)}w Zi t Zi ) 



(7) 



where t Zi is the time interval for redshift bin i and w Zi is the 
correction factor calculated from Equation (O. The full red- 
shift range is 0. 1 to 0.4 for the lowest range of Lx and 0. 1 
to 0.6 for the remainder. In the bins for each range of Lx, 
the evolution of Pj Ct seen in the right panel of Figure [7] is 
The cumulative jet power per cluster from jets more powerful overwhelmed by the large uncertainties, particularly from the 
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FIG. 7. — Left: Mean number of radio sources per cluster, per log Pi. 4 as a function of NVSS radio power, Pi. 4. Right: Cumulative jet power per cluster for 
radio jets more powerful than Pj ct . Details are given in §6\ Both functions are calculated for the four redshift ranges listed in the legend of the left panel. 



scatter in the relation of Equation (Q]i. Because of this, possi- 
ble differences between the serendipitous and all-sky surveys 
disc ussed i n iQare a minor issue. 

In MMN11, we concluded that (Pj C t) shows no significant 
dependence on Lx- However, the limited sample size pre- 
vented isolation of Lx from the redshift, because the most lu- 
minous clusters in the 400SD sample are at higher redshifts. 
Using the larger cluster sample here, we can break this degen- 
eracy and estimate the (-Pjct) for clusters with different X-ray 
luminosities over the redshift range 0.1 < z < 0.6 and the 
upper panel of Figure[S]shows that the increase of (Pj e t) with 
X- ray lumin osit y is not significant, consistent with the results 
in IMMN1 ll and iGiodini et alj d2010h . Since (P jet ) is simi- 
lar for all clusters, regardless of their X-ray luminosities, the 
energy input per particle from AGN is larger in less massive 
clusters. 

Using the L x - M 500 relation o flVikhlinin et al] (120091) and 
a gas mass fraction of 0.12, we can estimate the average gas 
mass within P500, (Mgas), for each luminosity range in the 
upper panel of Figure [8] Integrating the jet power over time 
(cf. Equation [7]), gives the mean energy injected into clusters 
by radio AGN. Therefore, the mean total energy per particle 
injected by the radio sources is 



t-! V jct/int [ -z,intM'"-p m\ 

Ejct = — — ' (8) 

where /1 = 0.59 is the mean molecular weight, m p is the 
proton mass and t Zy i nt = 'Yl,A !Si . Here, the integration time 
is limited to correspond to the redshift ranges for the sample 
bins. In principle, the energy injected by the radio sources 
sho uld be traced back to the time when BCGs formed (z ~ 2, 
e.g.. Ivan Dokkum &^ Franx 2001). Even with no AGN evolu- 
tion, extending the integration back to z = 2 boosts the energy 
injected per particle substantially (red arrows in Figure[8]l over 
the values for the redshift range of the sample (small dots). 
This estimate is conservative, because AGN are more active 
in the past (e.g jGalametz et al.l2009tlMartini et al.l2009l) and 
the clusters have assembled from smaller systems that may 
well have contained more than one BCG. To allow for the 



evolution of the Pj et Equation © can be generalized to 

where we assume that (M gas ) does not depend on the time. 
The evolution of the cumulative Pj C t (Figure [7j right) is mod- 
eled using a simple linear function, 

(Pjct)int ~ ^(Pjct) ~ A + Bt, (10) 

where the parameters (A, B) are fitted to $ z ( t )(Pj et ) from 
Equation © for z{t) = [0.15,0.3,0.5], and P jet = 3 x 
10 44 erg s~ 1 . This model raises the total energy input by an- 
other factor of 50% (blue arrows in Figured). 

We discuss the interpretation of Figure[8]in the next section, 
fJH] In short, the average AGN energy input to the clusters with 
luminosities of 0.3 < L^(10 44 ergs" 1 ) < 1.0 can reach 1.3 
- 2 keV/particle for ICM within P500, depending on the de- 
tails of AGN evolution. For the most massive clusters, with 
X-ray luminosities of Lx > 10 ergs -1 , the average AGN 
input energy is also significant, at 0.2 to 0.3 keV/particle. 
Note that the energy input of the single AGN outburst in the 
MS0735+7421 cluster is ~ 0.2 5keV per particle within the 
central 1 Mpc dGitti et al.ll2007l) . It is therefore plausible that 
a single, powerful AGN outburst can rival the integrated AGN 
energy input over time. 

8. DISCUSSION 

8.1. AGN Energy Input 

A few points need to be addressed regarding the results in 
Figure[8] First, the average jet powers are affected dispropor- 
tionately by the most powerful radio sources, which are the 
least likely to be background sources. As shown in the lower 
panel of Figure[8] the average AGN energy deposited per par- 
ticle for X-ray luminous clusters would be boosted by a factor 
of two assuming the scaling relation Equation (0} holds for 
the few sources with Pi. 4 > 10 42 ergs _1 . This would im- 
ply that a single, powerful radio AGN can be as important to 
heating atmospheres as the integrated pow er ou tput of radio 
sources over time. As discussed in section [3~2l Equation (Q~|i 
may be overestimating the Pj ct for the few powerful sources 
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FIG. 8. — Upper panel: Mean jet power (Pj e t) vs. X-ray luminosity. Gray 
circles correspond to the bins in z and Lx shown in FigurefT] For each range 
of Lx, darker shades of gray correspond to lower redshift bins. The red 
points with uncertainties are time averaged values, (Pj e t)int, for the same 
ranges of Lx as tne g mv points. The number of clusters hosting NVSS 
sources and the total numbers of clusters is given at the bottom of the panel 
for each range of Lx- Lower panel: Average AGN energy injected per par- 
ticle. The filled circles with error bars give the average energy per particle 
injected by radio jets, integrated over the redshift ranges marked in FigurefT] 
Jet powers are calculated from Equation {TJ, excluding the most powerful ra- 
dio sources with Pi. 4 > 10 42 erg s — 1 . The impact of the jets is extrapolated 
to z = 2 using two models: constant mean jet power (red arrows) and the 
linear evolution of Pj ct discussed in [|6](blue arrows). The effect of including 
Pj et for the most powerful radio sources is also shown, calculated using the 
saturated scaling relation as diamonds and the scaling relation of Equation (TJ 
as triangles. 

with Pi .4 > 10 42 erg s _1 in Figure[2] an issue that can only be 
resolved with deep X-ray imaging. The saturated scaling rela- 
tion gives average jet powers including all radio sources (dia- 
monds) that differ little from those obtained when the power- 
ful radio sources are excluded (circles). If the scaling relation 
does saturate, the small offsets between the spheres and dia- 
monds in the lower panel of Figure [8] show that excluding the 
most powerful sources makes little difference to the mean jet 
power estimates. 

On the other hand, our assumption that the cluster masses 
(M500) remain constant for the calculation of the mean en- 
ergy injected per particle fails to consider the hiera rchical as- 
sembly of clusters. Allowing for cluster growth, lHart et ail 



d201 ll) demonstrated that jet power from central radio AGN 
in clusters could increase by a factor of 10 per particle from 
z = 0.2 back to 1.1. At earlier times, AGN jets are more 
powerful and cluster progenitors are less massive. Therefore, 
the energy from radio jets accumulates more quickly per ICM 
particle in the building blocks of present day clusters. As less 
massive clusters assemble to form more massive ones, most 
of the excess energy is preserved. By ignoring this effect, we 
have certainly underestimated the total AGN energy accumu- 
lated in the ICM of massive clusters over their histories. One 
effect that might counteract this is that radio jets could break 
out of the atmosph eres of less massive ha los (e.g., the poorly 
confined sample of Cavagnolo et al. 2010), allowing some jet 
energy to escape. However, the huge mass of gas that will 
form the atmosphere of an incipient cluster is present from 
the outset and very few jets can escape from that. Jet energy 
escaping from a smaller atmosphere will be deposited in sur- 
rounding gas that is fated to collapse into the cluster. Unless 
the energy deposited by a jet is sufficient to unbind this gas 
from the cluster, it inevitably collapses into the cluster at some 
later time, carrying the excess energy along with it (apart from 
energy lost to radiation). 

Our calculations ignore energy lost from the ICM by X-ray 
radiation. If we are interested in the net energy gain from ra- 
dio jets, this must be taken into account. On the face of it, the 
upper panel of Figure [8] shows that the X-ray power radiated 
by the clusters in our sample exceeds the average power input 
from the central AGN for Lx > 3 x 10 44 ergs -1 . Our esti- 
mates of the mean jet power are conservative and only include 
central radio sources (projected within 250 kpc), when other 
radio sources may a ugment the total energy input significantly 
dStocke et al.ll2009l) . Nevertheless, it is likely that the ICM of 
the most luminous clusters suffers a net energy loss. It should 
be borne in mind that most of the X-ray power radiated by the 
great majority of clusters does not originate from a cooling 
core. Outside cooling cores it will take a very long time for 
the energy loss to have any noticeable impact. Central cooling 
tunes are available in the ACCEPT database (Cavagno lo et alj 
2009) for only 1 10 members of our sample, leaving us poorly 
placed to examine the net effect of jets on cooling core clus- 
ters. Notably, Figure[8]shows that, even with our conservative 
estimates for the energy input from central radio jets, clusters 
less luminous than Lx — 3 x 10 ergs -1 see a net energy 
gain. Thus, energy injected by central radio AGN accumu- 
lates in lower mass clusters, so that the integrated energy gain 
shown in the lower panel of the figure is mostly retained in 
these systems and has a significant impact on the ICM. 

As discussed in |j5j our sample shows increases in the frac- 
tion of clusters with central radio AGN for increases in both 
the X-ray luminosity and the redshift (Figure [5]). It is, there- 
fore, surprising that we do not see a more pronounced increase 
in the mean jet power with X-ray luminosity in the upper 
panel of Figure[8] The primary cause of this is the large scatter 
introduced by using Equation (Q]) to convert radio powers to 
jet powers. There is good reason to believe that mean jet pow- 
ers do increase with cluster luminosities (masses). However, 
the modest increase is buried by the scatter in the Pi .4 — Pj 0t 
relation. It is clearly desirable to find a more accurate way to 
estimate jet powers. 

8.2. Supermassive Black Hole Growth 

The integrated power output from radio- AGN at the cen- 
ters of clusters over the past ~ 10 Gyr implies substantial 
supermassive black hole growth. We have estimated the ac- 
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creted mass required to fuel AGN from the integrated AGN 
power output over time. We assume a conversion efficiency 
between accreted mass and mechanical jet power of r\ = 0.1, 
where Pj 0t = r]Mc 2 , and we ignore radiation loses. Integrat- 
ing the AGN mechanical energies shown in the bottom panel 
of Figure |8] over 5.7 Gyr (z — 0.6) gives an average accreted 
mass of 2 — 5 x 10 8 M Q per supermassive black hole. Ex- 
trapolating back to z = 2.0 over a look-back time of about 
10.5 Gyr, and assuming the modestly rising AGN power dis- 
cussed earlier, implies an average increase of 6 — 14 x 10 8 M Q 
per supermassive black hole. Note that, in hierarchically as- 
sembling clusters, this mass may be distributed among several 
black holes. These values are comparable to the black hole 
mass es of BCGs infer red from black hole scaling relations 
(e.g. lLauer et ai]l2007l) . which are thought to have been im- 
printed during the quasar era. Our result implies that normal 
AGN maintained over time by hot atmospheres may be as im- 
portant to supermassive black hole growth in BCGs as earlier 
and, presumabl y, much more rapid fo rmation processes (see 
the review in lMerloni & Heiri3l2012l) . It is conceivable that 
normal radio-AGN activity may give rise to black hole masses 
in exce ss of the mass exp ected from the Mbh-c relation for 
BCGs ( lLauer et "ail 12007b. Furthermore, 77 may lie w ell be- 
low 0. 1 dMerloni & Heinzll2008t iGaspari et al1l2012l) . If so, 
the inferred black hole growth rates may be even larger, lead- 
ing to the p ossibility of growing ultramassive black holes in 
BCGs (e.g., iMcNamara et al.ll2009t iHlavacek-Larrondo et alJ 
I2012bl) . 

9. SUMMARY 

We have combined eight surveys of X-ray clusters to com- 
pile a composite sample with 1032 clusters located in the 
area covered by the NVSS. For each NVSS radio source pro- 
jected within 250 kpc of a cluster center, we have estimated 
the mechanical power of i ts radio jet using the scaling re- 
lation, Equation (Q~|), from Cavag nolo et alJ (120101) . The jet 
power is weakly correlated with the X-ray luminosity of a 
hosting cluster, but the most powerful radio sources, with 
Pi. 4 > 1.4 x 10 42 ergs -1 , are all located in massive, cooling 
core clusters. The correlation is stronger if only the strong 
cooling core clusters with t coo \ < 1 Gyr are considered. We 
have also examined the distribution of radio source powers in 
cooling and non-cooling c ore clusters, using valu es of t coo \ 
from the ACCEPT project dCavagnolo et al.lf2009). Based on 
the modest number of our sample clusters in the ACCEPT 
database, radio sources in non-cooling core clusters are, in 
general, as powerful as those in cooling core clusters, except 
that the most powerful sources mostly appear in cooling cores. 

We have examined both the average radio power of clusters 
and the fraction of clusters with radio sources. The cluster 
sample is large enough to separate the dependence of the radio 
source fraction on redshift and cluster X-ray luminosity and 
we find that it increases moderately with both. The average 
power is also larger in more massive clusters and at higher 
redshifts. 

Finally, we have calculated the average AGN jet power 
using the scaling relation in Equation ((TJ (Cavagnolo et al. 
2010). This overestimates Pj e t for the few powerful radio 
sources in our sample with Pi. 4 > 10 42 ergs -1 , so that the 
average jet power would be dominated by these extremely 
powerful sources. Two approaches were used to solve this 
problem. In the first approach, the most powerful sources 
were simply excluded from the calculations, giving a lower 
limit on the average jet power. In the second approach, Pj 0t 



for the powerful radio sources was determined using a satu- 
rated version of the scaling relation, with the saturation level 
set empirically and saturated jet power based on the cavity 
powers of the five most powerful sources in our sample. The 
average jet powers determined using these two approaches are 
similar. In the upper panel of Figure[8] the average jet power 
is plotted against the cluster X-ray luminosity. Although the 
average jet power for the most luminous clusters is higher than 
for less luminous clusters, large uncertainties in their estima- 
tion make the differences insignificant. In general, the average 
jet power exceeds 3 x 10 44 erg s -1 even in the least luminous 
clusters, with 0.3 < Lx < 1 x 10 44 ergs -1 . Thus, the aver- 
age jet power exceeds the radiation output of the least massive 
sample clusters by an order of magnitude. 

The average jet power was integrated to redshift z — 2.0 
using two simple evolutionary models for the radio sources. 
For the first model, the radio power was taken to be con- 
stant and then the average AGN energy injected by jets ex- 
ceeds 1 keV per particle in the least luminous clusters, with 
0.3 < L x < 1 x 10 44 crgs -1 , and ~ 0.2keV per particle in 
the most luminous clusters, with Lx > 10 45 ergs -1 . Here, 
the number of gas particles was calculated using the gas mass 
within P500, determined from the L\— -M500 relation. For the 
second model, the average jet power was taken to be a linear 
function of the cosmic time and then, integrating to z = 2.0, 
the AGN energy input amounts to ~ 2keV per particle for 
clusters with 0.3 < Lx < 1 x 10 44 erg s -1 and ~ 0.3keVper 
particle for clusters with Lx > 10 45 ergs -1 . 

Existing X-ray data for our sample are inadequate to distin- 
guish the energy radiated by gas that would be significantly 
affected by radiative cooling. However, the total radiation out- 
put of the less massive clusters is small compared to the en- 
ergy input from AGN. If the energy injected by AGN is stored 
in these systems, we estimate that the AGN energy injected 
since z = 2.0 is significant for preheating of clusters. If so, 
rather than preheating, the effect of the AGN would be better 
described as "continual heating." In carrying out these calcu- 
lations, we have ignored the hierarchical assembly of clusters 
by assuming that the cluster masses are fixed. Since massive 
clusters assembled from less massive clusters, where the jet 
power per particle is larger, we expect that our estimates of the 
total AGN energy accumulated in massive clusters are low. 
We conclude that continual AGN energy input in the "radio 
mode" could well provide > 1 keV per particle in less mas- 
sive clusters, which approaches the excess energy required to 
account for observed departures from the s elf-similar scalin g 
relations that would be expected otherwise ( Wu et al. 2000). 

Lastly, we have estimated the mass that was accreted by 
supermassive black holes in BCGs to fuel their radio AGN 
and power their jets. Assuming that the jet power is related to 
the accretion rate by Pj e t = O.lMc 2 , for a typical BCG in our 
sample, the nuclear black hole would have grown by about 
10 9 M Q since z = 2. T his is compa r able to black hole masses 
for BCGs estimated by lLauer et al.1 (120071) . implying that that 
the fueling of radio AGN at the centers of hot atmospheres 
may be as significant as the earlier quasar era for the growth 
of supermassive black holes in BCGs. 
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